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We investigate the potential to observe double parton scattering at the Large Hadron Collider in pp — > WbbX — ¥ 
ivbbX at 7 TeV. Our analysis tests the efficacy of several kinematic variables in isolating the double parton 
process of interest from the single parton process and relevant backgrounds for the first 10 fb — 1 of integrated 
luminosity. These variables are constructed to expose the independent nature of the two subprocesses in double 
parton scattering, pp — ¥ IvX and pp — > bbX. We use next-to-leading order perturbative predictions for the 
' double parton and single parton scattering components of Wbb and for the pertinent backgrounds. The next- 

to-leading order contributions are important for a proper description of some of the observables we compute. 
We find that the double parton process can be identified and measured with significance S/\/"B ~ 10, provided 
i the double parton scattering effective cross section a c g ~ 12 mb. 

00 

1. Introduction 

i-C ' The successful operation of the Large Hadron Collider (LHC) and its detectors opens a new era in particle 
physics. The higher energies and larger luminosities at the LHC make it possible to explore new physics scenarios 
p> and to investigate unexplored aspects of established theories such as quantum chromodynamics (QCD). 

The standard picture of hadron-hadron collisions is shown on the left side of Fig. [1] One parton from each 
proton partakes in the hard scattering to produce the final state. The probability density for finding parton i in 
a proton with momentum fraction Xi and at the factorization scale fi is parameterized by the parton distribution 
function (PDF) fp(xi,fi). In this single parton scattering (SPS) scenario, the differential hadronic cross section 
neatly factors into: 

^ : , 

da%p S = Y J fp(xi,fi)f^(x' 1 ,fj.)da ij (x 1 ,x' 1 ,fi)dx 1 dx' 1 . (1) 

Qv^ . The "short-distance" partonic cross section d&ij is computed in perturbation theory, whereas the PDFs are 
f^i ' non-perturbative objects and must be extracted from experiment. 

The full description of hadronic collisions involves other elements including initial- and final-state soft ra- 
diation, underlying events, and multi-parton interactions. Double parton scattering (DPS) describes the case 
in which two short-distance subprocesses occur in a given hadronic interaction, with two initial partons being 
active from each of the incident protons. The general picture of DPS is shown on the right side of Fig. [1] Some 
evidence for DPS has been observed at the CERN Intersecting Storage Rings Q, the CERN SPS @], and more 
recently, at the Fermilab Tevatron 0, Q ■ 

Under the assumption of weak dynamic and kinematic correlations between the two hard scattering subpro- 
cesses, the general approach in these studies is to assume the differential hadronic cross section takes a factored 
form in analogy to Eq. (TTJ): 

do-p p PS = J2 J H p k (x 1 ,X2,HA,HB)H p il (x' 1 ,x' 2 ,fi A ,HB) 

xda 2J (xi,x' 1 , ^A)daki(x2,x' 2 , HB)dxidx2dx' 1 dx' 2 , (2) 

where m is a symmetry factor which is equal to 1 (2) if the two hard-scattering subprocesses are identical 
(non-identical). The joint probabilities Hp> k (xi,X2, fJ>A, /As) can be approximated as the product of two single 
PDFs: 

H l p ' k (x 1 ,X2,HA,^B) = fp{x 1 ,fl A )fp(x2,VB) ■ (3) 

Given that one hard-scattering has taken place, the parameter er c g measures the size of the partonic core. 
Typical values in phenomenological studies focus on the 10-12 mb region, consistent with measurements from 
the Tevatron collider 0,|||. 
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Figure 1: Schematic depiction of single parton scattering (left) and double parton scattering (right). Labels are as in 
the text. 

In an earlier study @, we investigated the DPS and SPS contributions at the LHC to the four-parton final 
state pp — > bbjjX in which a bb system is produced along with two jets j. We showed that there are characteristic 
regions of phase space in which the DPS events are expected to concentrate, and we developed a methodology 
to measure the effective size of DPS. Precise measurements of DPS at the LHC will provide insight into parton 
correlations, non-perturbative dynamics in hadron-hadron collisions, the structure of the proton, and parton 
distribution functions. 

In this talk, we examine another final state besides bbjj which may be a good candidate to observe DPS, 
namely the production of a W boson in association with a pair of bottom quark jets. The cross sections for 
both subprocesses, W and bb, are individually large, and the charged lepton from the W decay (along with 
the bottom quarks in the final state) provides a relatively clean signal to tag on. Our purpose is to establish 
whether double parton scattering can be observed as a discernible physics process in Wbb production at LHC 
energies. We remark that the Wbb final state is a significant background for the production of a Higgs boson 
H in the HW ± mode, where the Higgs boson decays into a pair of bottom quarks [6j, and that it can be a 
background in channels where new physics may arise such as in single top quark production [7j. Once DPS 
production of Wbb is observed, it would be interesting to assess its potential significance as a background in 
such searches. 

The rest of the talk is structured as follows. In Section^ we outline our procedure for computing the DPS and 
SPS contributions to Wbb at the LHC, and we discuss and evaluate backgrounds to the same final state. Section 
131 is devoted to the role of the large tt background. In Section 0] we focus on discrimination between the DPS 
and SPS contributions to Wbb. We study various single variable and two-dimensional kinematic distributions 
to bring out the DPS contribution more cleanly. Section [5] contains our summary. 

2. DPS and SPS contributions to Wbb production at the LHC 

We begin with the premise that there are DPS and SPS components of the same Wbb final state. Our aim 
is to try to pick out the DPS component and to study its distinct properties. Here we outline our method 
for computing event rates for Wbb from DPS and SPS as well as backgrounds at the LHC. We perform all 
calculations at a center-of-mass energy of y/s — 7 TeV. Event rates are quoted for 10 fb _1 of integrated 
luminosity. For the DPS case, Wbb production is computed using Eq. ([2]) where it is assumed that one hard 
scattering produces the W boson via the Drell-Yan mechanism, while the other scattering produces the bb 
system. Schematically, we can represent the partonic DPS process as: 

(ij W /± ) ® (kl -> bb) . (4) 

The individual SPS processes which make up the DPS process in Eq. (|4| are generated using the POWHEG 
BOX event generator |8l-ll0l] which includes next-to-leading order (NLO) QCD corrections for both, plus shower 
emission. The <g> symbol denotes the combination of one event from each of the W^ and bb final states. All 
events are produced using CT10 NLO PDFs [Hi- 
Ill the SPS production of Wbb, one hard-scattering produces the complete final state. The events from this 
process are also generated using the POWHEG BOX [l2j which implements the NLO calculation of Ref. fl3l |. 

Many standard model processes imitate the Wbb — ¥ bblv final state. In particular, we consider contributions 
from the following final states: 
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• Top quark pair production it where either (i) both t's decay semi-leptonically (denoted by tgte), and 
one of the charged leptons is missed or (ii) where one t decays semi-leptonically while the other decays 
hadronically (denoted by tgth) and only two jets are observed. 

• Single top quark production (tb, tb, tj and tj modes) where t —> W + b (t — > W~~b) 

• Wjj, where both light jets are mistagged as a b jets 

• Wbj where the light jet is mistagged as a & jet 

Other processes were found negligible after cuts. We compute the Wjj DPS contribution using the same method 
discussed above for Wbb, using the additional POWHEG code [l4j]. The top pair [l(| and single top SPS 
processes are also generated using the POWHEG BOX. Other SPS processes are generated using MadEvent [l6| 
or ALPGEN dT]; Wjj SPS is reweighted using a K-factor obtained with MCFM [Hj]. 

In order to avoid soft and collinear divergences in the processes that we compute using MadGraph and 
ALPGEN, we apply a minimal set of generator-level cuts: 

p Tj >15GeV, M<4.8, \ m \ < 2.5 , p T ,b > 15GeV , \r) b \ < 2.5 , AR mm > 0.4 , AR jm >0A, (5) 

where r\ is the pseudorapidity and Ai?^ = \J (rji — rjk) 2 + (</>/ — 4>k) 2 is the separation in the azimuthal- 
pseudorapidity plane between the two objects I and k. 

2.1. Simulation 

To identify the Wbb final state and reduce backgrounds, we begin with simple identification cuts on the 
generated event samples. First, we consider only muonic decays of the W boson (W — > fiv). We limit the 
hadronic activity in our events to include exactly two hard jets, both of which must be identified as bottom 
quark jets. Finally, all events (DPS and SPS Wbb as well as backgrounds) are required to pass the following 
acceptance cuts: 

PT.b > 20 GeV , \r) b \ < 2.5 ; 20 GeV < p T ,^ < 50 GeV , \n M \ < 2.1 ; E T > 20 GeV; AR bS > 0.4 , AR hfl > 0.4 . (6) 

We also apply detector resolution smearing and mistags; details can be found in our paper [l9| . 

Table 12.11 shows the number of events from the Wbb final state (DPS and SPS) and the backgrounds both 
before and after the acceptance cuts are applied. In these results and those that follow, we sum the W + and 
W~ events. In evaluating the DPS processes, we assume a value a c ff ~ 12 mb for the effective cross section, 
but the goal is to motivate an empirical determination of its value at LHC energies. The results in Table 12.11 
make it apparent that Wbb production from SPS and the top quark pair background are the most formidable 
obstacles in extracting a DPS signal. 



Table I: Numbers of events before and after the various cuts are applied for 10 fb 1 of data. 



Process 


Generator-level Cuts 


Acceptance Cuts 


p T < 45 GeV 


S' PT < 0.2 


W ± bb (DPS) 


10000 


247 


231 


173 


W ± bb (SPS) 


44000 


1142 


569 


114 


tt 


225000 


1428 


290 


13 


W ± jj (DPS) 


476000 


43.5 


37.7 


27.3 


W ± jj (SPS) 


20300000 


101 


55.7 


19.6 


Single top 


20000 


492 


168 


15 


W ± bj 


153000 


152 


53.1 


8.2 



3. tt Background Rejection 

We examine three possibilities to reduce the ti background: a cut to restrict fir from above, rejection of 
events in which a top quark mass can be reconstructed, and a cut to restrict the transverse momentum of the 
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leading jet. In the end, an upper cut on ,Et in the event appears to offer the best advantage. Indeed, one would 
expect that ,Et in Wbb events would be smaller than J£y in ti events. Top quark decays give rise to boosted 
W 8 which, after decay, should result in larger values of missing Ex compared to the Wbb process. The ,Et 
distribution is shown in Fig.[5]for the DPS component of Wbb, the SPS component of Wbb, and all backgrounds 
(left). On the right, we show the DPS component of Wbb and the ti background alone. The plot on the right 
shows that the DPS signal is produced in the region of relatively small f^x and the ti background has a harder 
spectrum in JCy. One way to suppress the ti background while leaving the DPS signal unaffected is to impose 
a maximum J?t cut in the 40-60 GeV range. In the analysis that follows, we include a maximum J£t cut of 45 
GeV in addition to the acceptance cuts outlined above. The effects of the maximum y^x cut are shown in the 
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Figure 2: The event rate as a function of Fjt for DPS and SPS. On the left, all backgrounds are included while the plot 
on the right compares the DPS events to those from ti alone. 



fourth column of Table [2TT1 This cut eliminates about 80% of the ti background that remained after the initial 
acceptance cuts. 

Backgrounds from events that contain a real top quark, such as the ti and single top events, might be reduced 
if one could reconstruct a top quark mass distribution Mm v from the final state objects (bottom quarks, charged 
leptons and neutrinos), and then eliminate events in which the reconstructed mass falls in a narrow window 
centered on the known top quark mass. We attempted this by asserting the reconstruction of the W mass to 
determine the longitudinal component of the neutrino momentum, however the ti background was not sufficiently 
peaked due to several issues. These are discussed in |19j . 

We also considered an upper cut on the px of the leading object, since this would tend to be harder for jets 
coming from a top quark. While the backgrounds did in fact have a larger px tail, they are also peaked on the 
low end where DPS events are located, and so we abandoned this observable as a possible discriminator. 

We conclude this section with the statement that of the three possibilities to reduce the ti background we 
considered, a cut to restrict fix from above appears to offer the best advantage and is the only added cut we 
impose in addition to the acceptance cuts specified above. 



4. Discrimination between DPS and SPS contributions to Wbb 



To separate the DPS events from those of SPS origin, we find it convenient to employ quantities which take 
into account information from the entire final state. One useful observable is S' pT , defined as [|| : 

9 , j_ If \PT(hM\ y , ( m*,") i y m 

PT V2VVMM + MW \\pt(1)\ + \pt(v)\J ' ( > 

Here px{b\, b%) is the vector sum of the transverse momenta of the two b jets, and px(£, v) is the vector sum of 
the transverse momentum of the lcpton and missing energy in the final state. 

In DPS production, the bottom quarks are produced roughly back-to-back such that the vector sum of their 
transverse momenta tends to vanish. Likewise, the vector sum of the lepton and neutrino momenta tends to 
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Figure 3: The event rate for Wbb production from DPS (left) and SPS (right) as a function of S' pri 



be small (with corrections from the boosted W ). Thus, the S' PT distribution for the DPS process exhibits an 
enhancement at low S' pT , as shown in Fig. [31 The peak does not occur at exactly S' pT — owing to NLO real 
radiation that alters the back-to-back nature of the bb and Iv systems. On the other hand, SPS production of 
Wbb final states does not favor back-to-back configurations, and it exhibits a peak near S' pT = 1. This feature 
is linked to the fact that many bb pairs are produced from gluon splitting Q. 

The clean separation in S' pT between the DPS and SPS Wbb processes exhibited in Fig.|3Jis obscured once the 
tt background is included (e.g., see the left side of Fig. [4]) . Figure [4] illustrates the effectiveness of the maximum 
]&t cut in reducing the ti background in the S' pT distribution. After the cut, a sharp peak is evident in the 
region of small S' pT where DPS events are expected to reside. 




Figure 4: The S' PT distribution for DPS and SPS production of Wbb including all relevant backgrounds. On the left, 
only the minimal acceptance cuts are imposed, while, on the right, an additional maximum cut is imposed (J5t < 45 
GeV). 



The plot on the right side of Fig.|U shows that extraction of a relatively clean DPS sample can be accomplished 
by imposing a maximum S' pT cut. The last column of Table HOI shows that a cut S' pT < 0.2 reduces the SPS 
Wbb rate while leaving the DPS signal relatively unaffected. In the end, the major background arises from DPS 
Wjj, as is expected since this process inhabits the same kinematic regions as the DPS Wbb signal. Despite this 
background, we find a statistical significance for the presence of DPS Wbb of Sj\f~B = 173/VT97 = 12.3. 
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4.1. Further discrimination 
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Figure 5: The event rate as a function of the angle between the normals of the two planes defined by the bb and iv 
systems (left) , and the azimuthal angle between the total momentum vectors of the bb and iv systems (right) . 



Observables which take into account the angular distribution of events are also useful in the search for DPS. 
In the left plot of Figure O we show the event rates for DPS Wbb and the backgrounds as a function of the 
angle between the normals to the two planes defined by the bb and Iv systems. These planes are defined in the 
partonic center-of-mass frame and are specified by the three-momenta of the outgoing jets or leptons. The angle 
between the two planes defined by the bb and iv systems is cos A9 b 5 £l/ = 713(61, 62) • n^{l, v), where h^(i,j) is 
the unit three- vector normal to the plane defined by the i — j system. We require full event reconstruction using 
the on-shell VF-boson mass relations to determine nz(l, v), and use both solutions to the quadratic equation. We 
see that the distribution of the DPS events is rather flat, aside from the cut-induced suppressions at Q b i tv ~ 
and ~ 7r, whereas the SPS events show a strong correlation, with a distribution that peaks near A0 b 5 tv ~ ?. 

In the right plot, we show the event rates as a function of the azimuthal angle between the total momentum 
vectors of the 66 and Iv systems. Since this azimuthal angle is defined in the transverse plane, it requires only 
,Et, and no neutrino longitudinal momentum reconstruction. The shape of the DPS distribution is flat while 
the SPS distribution shows a strong correlation, with a preference for values toward w. 

We also considered the azimuthal angle between the charged lepton and the total momentum vector of the 
66 system, since this was less sensitive to detector effects and higher-order corrections than using the neutrino- 
lepton system. However, this is far less efficient a discriminator, and for brevity, is omitted. 

In both plots of Fig. [SJ it is clear that DPS and SPS exhibit different behaviors as a function of angular 
observables. However, the dominance of SPS Wbb and backgrounds over DPS Wbb for the full range of these 
observables makes it impossible to extract a DPS Wbb signal from these distributions by themselves. 



4.2. Two-dimensional distributions 

Despite the dominance of the Wbb SPS contribution and the backgrounds over the DPS Wbb contribution, 
the angular distributions can still be extremely useful when used in conjunction with other observables. Two- 
dimensional distributions of one variable against another show distinct regions of DPS dominance (or SPS and 
background dominance). In Fig. [51 we construct one such scatter plot. We show S' pT versus the angle between 
the normals of the two planes defined by the 66 and Iv systems (AO^ lv ). We see that the DPS events reside 
predominantly in the lower half of the plane (small S' PT ) and are distributed evenly in the angular variable. The 
separation between DPS Wbb and the SPS component is not as pronounced in the S' pT — AS^ tv plane as we 
saw in our earlier study of bbjj pj. In the Wbb case, the background events are more evenly distributed over the 
full plane, to some extent resulting from inclusion of both solutions for the neutrino's longitudinal momentum 
in the W^ decay. 

In Fig. [7J the SPS Wbb and background events in the S' pT — A<j> b i e ^ T plane show a strong preference for 
the upper right-hand corner of the plane. This two-dimensional distribution indicates that cuts on the S' and 
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Figure 6: Two-dimensional distribution of events in the variables S' PT and A9jj ifr . The Wbb DPS events lie in the lower 
half of the plane, while the Wbb SPS and background events occupy the upper half. (DPS events are marked as 'x'; SPS 
+ backgrounds with a dot.) 

A0 b 5 i variables should permit extraction of an enriched sample of DPS Wbb events. To quantify the degree 
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Figure 7: The two-dimensional distribution of events in the variables S' PT and A<f>bb,8? T (left). The box denotes the 
boundary which gives the highest statistical significance. On the right, we show the significance dependence as a 
function of the corners of the box. 

of separation, we define a region in this plane that gives the highest statistical significance. Its boundary is 
denoted by the black box in the left panel of Fig. [JJ Restricting S' PT < 0.25 and A(j>bb,e~$ T < 37r/4, we find a 
statistical significance of S/V~B = 15.2. 

By employing distributions in both SL and Afibb^fyri we achieve a better significance than from S' pT alone. 
In the right plot of Fig. we show the dependence of the significance on the placement of the box. As long as 
the maximum value of A4>bb,£^ T 1S m the 7r/2 — 37r/4 range, a statistically significant extraction of DPS Wbb 
from the other events can be obtained, given our assumed effective cross section a c s = 12 mb and luminosity. 

We suggest experimental analyses of Wbb at the LHC in terms of the two-dimensional distributions presented 
in this section with the goal to establish whether a discernible DPS signal is found. The enriched DPS event 
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sample can be used for a direct measurement of the effective cross section a e ff. 



5. Conclusions 

In this talk, we investigate the possibility to observe double parton scattering at the early LHC in the 
pp — > WbbX — > IvbbX process. 

After identifying the most relevant background processes, we pinpoint a set of observables and cuts which 
would allow for the best separation between the DPS Wbb signal and the backgrounds (including the SPS Wbb 
process). To provide the most precise predictions possible, we generate the DPS Wbb signal event sample, 
the SPS Wbb sample, and the dominant background event samples at next-to-leading order in QCD. The 
main obstacles in the extraction of the DPS signal are the backgrounds from tt production and the SPS Wbb 
component. The most efficient way to suppress the tt background is with an upper cut on the missing transverse 
energy of the event, since top quark decays result in larger values of Et. 

To separate the DPS component of Wbb from the SPS component, we find it useful to employ observables 
which take into account information on the full final state rather than observables which involve one or two 
particles. Examples are the S' PT variable (defined in Eq. {2}) and the angle (AQ^ ^) between the two planes 
defined by the bb and Iv systems, respectively. By displaying the information from these two observables in 
two-dimensional distributions, we show in Sec. 14.21 that it is possible to identify distinct regions in phase space 
where the DPS events reside. Utilizing cuts on these observables that enhance the DPS Wbb sample, we find 
that the DPS signal can be observed with a statistical significance in the range S/^/B ~ 12 — 15. 

Once DPS production of Wbb is observed, it will be interesting to assess its potential significance as a 
background in searches for other physics, such as Higgs boson production in association with a W boson (where 
the Higgs boson decays as H — > bb), and precise studies of single top quark production where new physics could 
contribute to the Wtb vertex. 
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